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ABSTRACT 

In this report we develop an approximation scheme for the identification 
of hybrid systems describing the transverse vibrations of flexible beams with 
attached tip bodies. In particular, problems involving the estimation of 
functional parameters (spatially varying stiffness and/or linear mass density, 
temporally and/or spatially varying loads, etc.) are considered. The 
identification problem is formulated as a least squares fit to data subject to 
the coupled system of partial and ordinary differential equations describing 
the transverse displacement of the beam and the motion of the tip bodies 
respectively. A cubic spline-based Galerkin method applied to the state 
equations in weak form and the discretization of the admissible parameter 
space yield a sequence of approxirating finite dimensional identification 
problems. We demonstrate that each of the approximating problems admits a 
solution and that from the resulting sequence of optimal solutions a 
convergent subsequence can be extracted, the limit of which is a solution to 
the original identification problem. The approximating identification 
problems can be solved using standard techniques and readily available 
software. Numerical results for a variety of examples are provided. 

To be submitted for publication in the Journal of Mathematical Analysis and Applications. 
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SECTION 1 


INTRODUCTION 


In this report we develop an approximation scheme for the identifi- 
cation of systems describing the transverse vibration of flexible beams 
with attached tip bodies. The equations of motion for this type of prob- 
lem generally take the form of a hybrid system of coupled partial (gov- 
erning the vibration of the beam) and ordinary (describing the motion of 
the tip bodies) differential equations with appropriate geometric bound- 
ary conditions and initial data. The resulting identification problems* 
therefore, tend to have infinite dimensional constraints. Moreover, if 
the parameters to be identified are functional (spatially and/or tempor- 
arily varying), the admissible parameter space is a function space and as 
such is infinite dimensional as well. The solution of the resulting 
constrained optimization problem, therefore, necessitates the use of some 
form of finite dimensional approximation. 

The scheme we develop here is based upon the reformulation of the 
equations of motion in weak form. A cubic spline based Rayleigh-Ritz- 
Galerkin method is used to define a finite dimensional approximation to 
the state equations. Using finite dimensional subspaces to discretize 
the admissible parameter space, we obtain a doubly indexed sequence of 
approximating identification problems. Using standard variational argu- 
ments, we derive a convergence result for the state approximations. We 
show next that each of the approximating identification problems admits a 



solution and that from the resulting sequence of optimal parameter values 
a convergent subsequence can be extracted whose limit is a solution to 
the original infinite dimensional identification problem. In addition to 
convergence results, we present numerical results which demonstrate the 
feasibility of our method. 

The approach described in this report represents a significant 
improvement over the method developed in [22]. Indeed, we have developed 
a scheme which is computationally simpler and, by relaxing the necessary 
hypotheses on the admissible parameter space, is applicable to a wider 
class of problems. Our results are similar in spirit to those presented 
in [7] in the context of parabolic systems, in [8] for hyperbolic systems 
and in a forthcoming paper by Banks and Crowley [5] for beam equations 
with standard boundary conditions (e.g., clamped, simply supported, 
cantilevered, etc.). Other work regarding approximation methods for 
inverse problems in elasticity can be found in [2], [3], [4], [10] and 
[14]. 

We simplify our presentation by only considering a cantilevered 
beam with an attached tip (point) mass. As is discussed in Section 4, 
however, our general approach is applicable to a broad class of beam-tip 
body vibration problems. In Section 2 we derive the weak form of the 
equations of motion, define weak and strong solutions, and formulate the 
identification problem. In Section 3 we define the approximation scheme 
and discuss convergence. Numerical results are presented in Section 4. 

We employ standard notation throughout. The Sobolev spaces of real 
valued functions on the interval [a,b] whose kth derivatives are L 2 

are denoted by Hk(a,b), k = 0,1,2 The corresponding Sobolev inner 

products and their induced norms are denoted by <• ,.> k and | • | ^ 
respectively. For Z a normed linear space with norm J • j 2 and 

f : [0,T] Z, we say that feL 2 ([0,T],Z) if /J jf (t ) j ^ dt < «. 
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Similarly, f will be said to be an element in C^([0,T],Z) if the map 
t -> f(t) from [0,T] into Z is z times continuously differentiable on 
(0,T). Finally, for a function of one or more real variables the symbol 
D 0 f (D 0 f) will be used to denote the l s ^ (kth) derivative of f with 
respect to the independent variable 0. If f is a function of a single 
variable only, the subscript may be deleted. On occasion the short-hand 
notation D e f(eo) or Df(eo) will be used in place of D e f J Qo or Df j 0 q to 
denote the derivative of f evaluated at 0 q. 
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SECTION 2 


THE IDENTIFICATION PROBLEM 


We consider a long slender beam of length i with spatially varying 
stiffness El and linear mass density p which is clamped at one end and 
free at the other with an attached tip mass of magnitude m (see 
Figure 2.1). 



Using the Euler-Bernoulli theory and elementary Newtonian mechanics, we 
obtain the partial differential equation 

p(x)D*u(t,x) + D^EI (x)D^u(t,x) 

D x a(t,x)D x u(t,x) + f(t,x) xe (0,£) 

and boundary condition at the free end 


te(0,t) 


( 2 . 1 ) 
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( 2 . 2 ) 


mD 2 u(t,£) - D x EI(£)D^u(t,£) = 

- cr(t,*)D x u(t f ji) + g(t) te(0,T) 

describing the transverse displacement of the beam and tip mass 
respectively where a is the internal tension (as a result of axial 
loading) f is a distributed lateral load applied to the beam and g is a 
force directed transversely which acts on the tip mass (see [9] and 
[26]). Rotational equilibrium at the free end yields 

D x u(t,£) = 0 te[0,T] , (2.3) 

while at the clamped end we have the usual geometric boundary conditions 
u(t,0) = 0 Du (t ,0) = 0 te[0,T]. (2.4) 

A 

The temporal boundary conditions, or initial conditions, are assumed to 
be of the form 

u(0,x) = <f,(x) D t u(0,x) = *(x) xe[0,£] . (2.5) 

We make the standing assumptions that m > 0, EI,peL (0,£) with EI,p > 0, 

oo 

ae ^([OJLH^O,*)), g e L 2 (0,T), f e L 2 ([0,T],H°(0,T)) , <{, e H 2 (0,£), and 

^eH°(0,Jt ) with i|)(£) specified in R. Define the Hilbert space 
H = R x H^(0,£) with inner product 

< (n»<i>), (c»ip) > H = n? + < >q. 
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We then rewrite Equations (2.1) - (2.5) as 


. 2 * 


M 0 D t U (t) + A o u(t) = + F 0 {t) te(0,T) 


Y 0 u(t) 


= 0 Yl u(t)| = 0 
x=0 1 x=0 


T20<t) |x=« 


= 0 te[0,T] 


u ( 0) = I D t u(0) = $ 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


where u(t) = (u(t,£),u(t,.))eH, F Q (t) = (g(t ) ,f (t, • )) , $ = (<(> ( Jt ) , cj>) , 

4> = ($(£)>$) and the operators M Q , A Q , B Q (t) and , i=0,l,2 on H are 
defined formally by 

Mq (n ,<p) 

Aq (n ,<(>) 

B Q (t)(n ,<j>) 

Y-j (n,<f>) 

respectively. 

There exist several ways in which the notion of a solution to the 
system (2.6) - (2.8) can be made precise. Of particular interest to us 
here are the ideas of a weak or variational solution and a strong 
solution. 


= (mn »p4> ) » 

= (- DEI(£)D 2 <f>U), D 2 EID 2 4») , 

= (- <r(t,£)D<j>(fc), DaD*), 

= D 1 (j, , i = 0,1,2 
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Define the Hilbert space {V, < • > v } by 


V = { (n><|>)eH : < t »eH 2 ( 0, Ji ) , <f,(0) = D#(0) = 0, n = ♦ (*)} , 
< (<f>U)»<J>) , (<|>U),<J>) > v = < D 2 <{> , D 2 ip >q. 


It is not difficult to show that V can be densely embedded in H. 
Choosing H as our pivot space, we have therefore that VchcV' where V' 
denotes the space of continuous linear functionals on V. Consider the 
second order initial value problem 


< MgD^u(t), 0 > H + a(u(t), e) = 
b(t)(u(t),0) + < F 0 (t), e > H te (0,T) , 0eV 


(2.9) 


u(0) = $ D t u(0) = 


( 2 . 10 ) 


where 0 = (0(t),e) and the bilinear forms a: VxV -»■ R and b(t): VxV -> R 
are given by 


a($,(|>) 


< E I D 2 cf> , D 2 f > 


0 


and 


b(t)(<J>,<^) “ ^ " oD<j), Dip >q 
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respectively. A solution u to (2.9) and (2.10) with u(t)eV is known as 
weak or variational solution to (2.6) - (2.8). Indeed, if the deriva- 
tives in (2.6) and (2.7) are taken in the distributional sense, A Q 
and Bg(t) become bounded linear operators from V into V' with 

< A Q $,$ > h = a($,$) 


and 


< B 0 (t)$,$ > H = b(t)($,$) 

where the H inner product is interpreted as its natural extension to the 
duality pairing between V and V' (see [1], [19], [24]). Since 
Fg(t)eHcV' we have therefore that the systems (2.6) - (2.8) and (2.9) 
and (2.10) are two representations for the same initial value problem in 

V'. 

Under the assumptions which we have made above, standard arguments 
(see [16], [17]) can be used to demonstrate the existence of a unique 
solution u to (2.9) and (2.10) with ueC([0,T],V), D ueC([0,T],H) and 

D^ e L 2 ([0,T],V). 

In order to characterize strong solutions we rewrite (2.6) - (2.8) 
as an equivalent abstract first order system and then rely upon the 
theory of semigroups and evolution operators. Let Z = V x H with inner 
product 


< (vj.hj), (v 2 ,h 2 ) > z = a(v 1 ,v 2 ) + < Mghj.hg > H 
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We assume that EIeH^(0,jt) and aeC*([0,T], H*(0,£)) and define the opera- 
A: Dom(A)c Z + 1 by 


Dom(A) = Doiti(Aq) x V 


A 


-"o 1a o 


I 

0 


where I is the identity on V, Mg and Ao are as they were defined 
above, and 

Dom(Ag) = {$ = (<f>U), <f>) eV: <|> e H 4 (0,£), D 2 <j>(a) = 0} . 

Similarly, define the operators B(t): Z -»• Z by 


B (t ) 


0 0 

M^Bgtt) 0 


and let A(t): Dom(A)cZ + Z be given by A(t) = A + B(t). Let F(t)eZ be 
defined by F(t) = (0, Mg^Fg(t)), ZqeZ by Zg= and consider the 

initial value problem 

D t z(t) = A(t)z(t) + F (t ) te (0,T) (2.11) 

z(0) = z Q . (2.12) 
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It is not difficult to argue that the operator A is densely defined and 
conservative. That is 


< Az,z > z =0 zeDom(A). (2.13) 


Moreover, we have 

Theorem 2.1 : The operator A:Dom(A)cz -► Z is skew self adjoint. 

Proof 

We first argue that -Ac A*. That is Dom(A)c Dom(A*) and for 
zeDom(A) , A*z = -Az. Let z 15 z 2 eDom(A) with z. = ). Then 

< Az^ ,z 2 >2 - ^ »<i>2 ) + ^ ~ Aq$^ » 4 > 2 

- < EID D <f> 2 >q - < EID ,D \\>2 >q 

- ^ ^ 1 *” ^ 0^2 ^ 

= - <z 1 ,Az 2 >2 

where we have used integration by parts, the definition of Dom(Ao) and 
the definition of V in performing the above computation. This, of 
course, implies that z 2 eDom(A*) and A*z 2 = -Az 2 . 

We next argue that Dom(A*)c Dom(A). Let weDom(A*) and y = A*w. 
Then for zeDom(A) 


< z,y > z = < z,A*w > z = < Az,w > z - 
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Recalling that z,w,y e Z, let z = (z^Zg), w = (w p w 2 ) and y = (y^). 
Then 

0 = <z,y> z - <Az,w> z 

= aCzj,^) + < M 0 z 2 ,y 2 > H - a(z 2 ,w j) + <A q z 1 ,w 2 > h 
= <EID 2 z p D 2 y 1 > 0 + mz 2 (£)y 2 + < P z 2 ,y 2 > Q - 

<EID 2 z 2 ,D 2 Wl > 0 - DEI(£)D 2 z 1 (£)w 2 + <D 2 EID 2 z lf w 2 > 0 (2.14) 

a a 1 a 2 a 2? 2 

where y 2 = (y 2 ,y 2 )eH and w 2 = (w 2 ,w 2 )eH. Let 0eH (0,£) be defined by 

D 2 e = py 2 , e(£) = 0 and De(£) = - my 2 . 

Then substituting into (2.14) and integrating by parts, we obtain 

0 = - DEI(£)D 2 z 1 (£)(y 1 (£) + w 2 ) + 

< D 2 EID 2 z 1 ,y 1 + w 2 > Q + < D 2 z 2 , e - EID^ > Q 

which implies 

(i) -DEI(£)D 2 Zl (£)( yi (£) + w*) + <D 2 EID 2 z 1 ,y 1 + w 2 > Q = 0 
and 

(ii ) < D 2 z 2 ,9 - EID 2 Wj > 0 = 0 . 
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Let (n,$) be an arbitrary element in H. Choose zieH 4 (0,£.) which 
satisfies 


2 2 
D tID Zj 


♦ . DEI(£)D 2 Zl (£) = - n , D^Zj (i) 


0, Dz^(O) = 0 


and 


z l(0) = 0. 

Then z\ = (zi(£),zi)eDom(Ao) and (i ) therefore implies that 

< (yjU) + w^yj +W 2 ), (n »<f> ) > H = ° (n,+ )eH 

from which it follows that yjU) = -wJ, and y 1 = -w ^ in H°(0,£). We 
have, therefore, that = -y^eV. 

Next, let $ be an arbitrary element in h 0(0,£). Choosing 

z ? = / / ■}>, (ii) implies that 
0 0 

< e - EID 2 Wl , 4 > >0 = 0 (J>eH°(0,jz ) 

and hence that 0 = EID^w-^. This in turn implies that WjeH^( 0,£) and 
2 1 

0 w^U) = Hjjy 0 ( £ ) = Since weZ, w^eV and therefore WjeDom(Ag), we 
conclude that w = (w^ ,Wg )eDom(AQ) x V = Dom(A) and the theorem is proven. 
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Since A is densely defined [12, Theorem 3, page 142] implies that A* 
is closed and, by Theorem 2.1 above, that A is closed as well. This fact 
together with (2.13) and Theorem 2.1 yield that A is maximal dissipative, 
and hence by [13, Theorem 4.2, page 84] that it is the infinitesimal 
generator of a Cq semigroup of contractions { S (t ) : t >. 0} on Z. It is 
in fact the case that (see [13], [28]) S(t) is defined for t < 0 and that 
{ S (t ) : - °° < t < forms a Cq group of unitary operators on Z. Since 

the operators B(t), 0 <_ t <_ T are bounded [20, Theorem 2.3, page 132] 
implies that {A(t)}t e [0,T] is a stable family of infinitesimal 
generators. Since Dom(A) is independent of t, [20, Theorem 4.8, page 
145] yields that the family {A(t)}t e [o,T] generates an evolution system 
{ U ( t , s ) : 0 _< s _< t _< T} on Z. 

Def i ne 


z(t) = (u(t),v(t)) 


(2.15) 


in Z by 


t 

z(t) = U(t,0)z n + / U(t,s)F(s)ds. (2.16) 

u 0 

The continuous function z given by (2.16) above is the unique mild 
solution to the initial valve problem (2.11), (2.12). If, in addition, 
zoeDom(A) (that is $eDom(Ao) and <j>eV) then z is a strong solution to 
(2.11), (2.12). Indeed z is differentiable almost everywhere on [0,T] 
D 2 Ze ^2^ E0»T]»Z) » satisfies (2.12) and satisfies (2.11) for almost 
every te[0,T] and is such that z(t)eDom(A) a.e. on [0,T]. 


13 



We shall call u given by (2.15) and (2.16) a strong solution to the 
initial value problem (2.6) - (2.8). The following result is easily 
obtained. 

Theorem 2.2 : Suppose EIeH^(0,t), aeC*([0,T], H*(0,&)), $eDom(Ag), and 

JjjeV. Then u given by (2.15) and (2.16) is the unique strong solution to 
the initial value problem (2.6) - (2.8). We have that u satisfies (2.8) 
and (2.6) and (2.7) a.e. on [0 ,T] . Moreover, u is twice differentiable 
in H and differentiable in V almost everywhere on [0,T] with 
D2ueL 2 ([0,T],H) and D t ti e L 2 ([0,T],V) . 

A 

It is also not difficult to show that if a strong solution u to 
(2.6) - (2.8) exists it coincides with the weak solution, and in which 
case, it is given either by (2.15) and (2.16) or as the solution to the 
initial value problem (2.9), (2.10). 

In formulating the identification problem, for ease of exposition, 
we consider a (reasonably broad) class of inverse problems which are of 
particular interest in the structural modeling of large flexible 
spacecraft and shuttle attached payloads (see [25], [27]). We assume 
that we are only interested in estimating the parameters itieR, El, 
peLco(0,t), and aeL2([0,T],H 1 (0,t)) where a = a(ao,m,p) with 
aoeL2(0,T). It is not difficult to further generalize the results 
which follow to allow for the identification of initial data, the 
external loads f and g, and more general forms of the internal tension a 
(see [7]). The motivation for choosing a to be a function of a time 
varying function ao, the magnitude of the tip mass m,and the linear 
mass density of the beam p will be made clear below. 
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Let 21 = R x L<»(0,j!, ) x L«(0,O x L2(0,T) with the usual product 
topology. Let Q be a subset of 21 which satisfies 

(HI) Q is compact. 

(H2) There exist constants m-j, M-j , i=l,2,3 such that 

0 < <_ m _< 

0 < m 2 < El < M 2 

0 < m^ < p < M^ 

for all q = (m, EI,p,ao)eQ. 

(H3) For all qeQ, a(q)eL2([0,T]),H 1 (0,j> )) with the mapping 
q a(q) continuous from Q into L2([0,T],H 1 (0,fc)) . 

We assume that we have been provided with displacement measurements 
{U(t. ,x. )} *_i for the beam at positions x.e[0,£], j = 1,2, ...v, at 

j —1 * • • • v 

times t-je[0,T], i = 1,2, ...y, which result from a known input distur- 
bance applied to the system in a known initial state and formulate the 
identification problem as a least squares fit to data: 

(ID) 

Find q = (m,EI ,p ,ao)eQ which minimizes 

J(q;u(q) ) = \ l I u”(t, ,x, ) - u (t. ,x . ,q) 1 2 (2.17) 

i=l j=l 1 1 J 1 J 1 

where u(t;q) = (u(t,Ji;q),u(t,» ;q)) is the solution to (2.9), (2.10) 
corresponding to qeQ. 

The infinite dimensionality of both the state, which is governed by 
the system (2.9), (2.10), and the admissible parameter space Q (being a 
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function space) necessitates the use of some form of approximation in 
solving Problem (ID). We develop and analyze one particular scheme in 
the next section. We note that the approximation theory to be developed 
below will also permit the formulation of the identification problem 
based upon criteria other than displacement; for example, velocity (see 
[ 22 ]). 
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SECTION 3 


AN APPROXIMATION SCHEME 


Our scheme is based upon the construction of a sequence of 
approximating identification problems, in each of which both the state 
constraint and the admissible parameter space are finite dimensional. We 
argue that each of the approximating problems admits a solution. The 
resulting sequence is shown to have a sub-sequential limit which is a 
solution to the original identification problem. The state approximation 
is constructed using a spline based Galerkin method. The admissible 
parameter space is discretized using splines as well. We begin by 
discussing the state approximation. 

N 

Working abstractly at first, let V c V be a finite dimensional 

N N 

subspace of H and let P denote the orthogonal projection of H onto V 

with respect to the H inner product. The Galerkin equations 

corresponding to (2.9), (2.10) are given by 

< M q d£u (t), e > H + a(u (t),e ) = 

(3.1) 

b(t)(u N (t),0 N ) + < F Q (t),e N > H te (0,T) , 0 N eV N 

u N (0) = P N $ D t u N (0) = P N $ (3.2) 

m N N N 

where u (t) = (u (t,£),u (t,»))eV . We define the following sequence of 
approximating identification problems. 
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(IDN) 

Find q = (m,EI,p ,a Q )eQ which minimizes J(q;u N (q)) where J is given 

by (2.17) and u N (t;q) = (u N (t ;q ) ,u N (t , . ; q)) is the solution to 
(3.1), (3.2) corresponding to qeQ. 

Let } j=-i de note the standard cubic B-splines on the interval 

[0,£] corresponding to the uniform partition a N = {0, ^ , ^,...£} (see 
[21]). Let denote the modified cubic B-splines which satisfy 

Bj (0) = DBj(O) = 0, j = 1,2, ..N+l. That is 


bJ = bJ| - 2bJ - 21^ 


B N = 

J J 


j = 2,3. ..N+l, 


*N N N N 

Let B. = (B.(t),B.) and let V c V be defined by 

J J J 


N -N N+l 

V = SPAN {B j} j=1 


The Galerkin equations (3.1), (3.1) take the form 


Mq*W^ (t ) + AqW N (t ) = Bg(t) W ^ (t) + Fg(t) 


(3.3) 


w"(0) = (W N )" 1 $ N w» = (W N )‘ A $ 


Nx-UN 


(3.4) 
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where 


Klij 

d N /o x d N,„x . / d n d n 
= mB . (t)B •(£,) + / pB.B. 

' J Q * J 


= /eid^dW 
0 1 J 


= -/a(t,.)DBM 
0 J 

[■>>]< 

= g (t )B^( JZ.) + J f (t,» )B? 
1 0 1 

[i N ii 

= 4> (5. ) + J ^ 

1 0 1 

[i N ], 

£ 

= tJj(^)B^U) + / ifiBV* 

1 0 1 


■ b M*> + 

-N, X 

and u (t) = 

N+1 N aN 

& "J (t) V 


Our convergence arguments are based upon the approximation 
properties of spline functions. Let So(a N ) = SPANiB^}^' 1 , and let 

SV) = SPAN{B^}^. For <f> a function defined on the interval [0,£] let 

J J A 

— N — 3 N 

I <}» denote that element in S (a ) which satisfies the interpolatory 
constraints ("i%)(^) = » j = 0,1,2...N, 0(7%) = 0<|>(j^-) j = 0 
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and N and let I^<(> denote that element in s3(a^) which satisfies the 
interpolator constraints (I N <f>)(-jjp) = j = 1,2. ..N, D(I N <{,)(£) = 

D<j> ( jz. ) . The interpolator spline will be well defined whenever <f> is 

well defined at the node points and D<f> at the end points. A similar 

N 

statement can be made for I <)>. 

We require the following two standard results concerning the 
approximation properties of interpolator splines (see [23]). 

Proposition 3.1 : For <|>eH 2 (0,£) 


|D k (A - ^)| 0 l cl k N- 2+k |D^| 0 k = 0.1 

where is independent of <}> and N. 

Proposition 3.2 : For <f>eH4(0,£) 

|D k (A - *)| 0 < cjV 4+k |D%| 0 k = 0,1,2 

2 

where is independent of and N. 

Lemma 3.1 

-N -N Na N N 

(1) Let <|> = ( 4 > ( £ ),<f>)eV and let cf> = P 4> = (<|> U),<j> ). Then 

N 2 ■'N A 

<f> (J) in H (0,£ ) and consequently (f> <(> in V. 

(2) P^ -► I strongly in H. 
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Proof 


0) 

|* N -+|ol |? N - *|h < |(A) -+|„ 

* | A -*| 0 . |A - *| 0 

<_ cJn" 2 |D 2 *| 0 » 0 as N » » 

A 

where (I N $) = ( (I N *)U), 1%). 

|D(4> N - <J>)| 0 < |D(/ - A)| 0 + |°(A - <t> ) | o 

< kjN J <j> N - 7\| 0 + cJn _1 |D 2 4)| 0 

< k 1 N|cj> N - *| Q + k lN |T\ - 4> J o + C^N -1 jD 2 <j,| 0 

< k 1 C^N“ 1 jD 2 <t , j Q + k^N -1 ^ 2 #^ + C 1 N" 1 |D 2 ({,| 0 

= (2^^ + c|)N -1 Jd 2 <{) J 0 -v 0 as N ->- « 

where we have made use of the Schmidt inequality (see [23]) in making the 
estimates above. 

Using the Schmidt inequality together with the first integral 
relation (see [23]), we obtain 
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|D 2 <J) N | Q < 2|D 2 (/ - 7%)| 2 + 2|D 2 A|o 

< 2k 2 N 4 |/ - A|o + 2 | D ^|o - 2 | D *(4» - A)|o 

< 2k 2 N 4 |/ - *| 2 + 2k 2 N 4 |T fi 4> - «^| 2 + 2|D%| 2 

< a 2 |D 2 «},| 2 

where a is independent of <}> and N. Let ipe { 0 : 0eH4(O,£),e (0) = De(0) = 
0} . Using arguments similar to those used above together with 
Proposition 3.2, it can be shown that 

|D 2 (/ - *)| 0 1 kN" 2 |D\| Q 

where k is independent of ip and N. Then 

j D 2 (4> N - 4> ) J 0 <_ jD 2 (<j) N - ) | 0 + JD 2 (i> N - t|»)j 0 + JD 2 (iJ; - ) J q 

<.(<* + 1)|D 2 U - <}> ) j q + kN 2 |D%Jq. 

Standard density arguments guarantee that tp can be chosen so as to make 
the first term arbitrarily small and therefore that 

|D 2 (<j, N - $)| Q + 0 as N * » 

which proves statement (1) of the theorem. 


22 



Statement (2) follows from Statement (1) and the fact that V is 
dense in H and the PN being orthogonal projections are uniformly 
bounded. 

Theorem 3.1 : Let {q^}cQ with + q as N ®. Let u denote the 

solution to (2.9), (2.10) corresponding to q and let u N denote the 
solution to (3.1), (3.2) corresponding to q^. Suppose further that u is 
a strong solution (see Theorem 2.2). Then u N -»• u in V and D t u N ->• D t u in 
H as N -*■ °° for each te[0,T]. 

Proof 

Our argument is similar in spirit to the one used in [11] to 

demonstrate the convergence of a Galerkin method for the integration of a 

class of hyperbolic systems. We adopt the convention that the 

superscript N on a form or operator indicates that it be computed with 
N N N N N 

respect to q = (m ,EI ,p ,3 q) while no superscript indicatesthat it be 
computed with respect to q = (m,EI,p ,a^) . 

Applying the triangle inequality, we have 

1 A N A | . I *N d N' , | . I rvN* A | /0 c v 

| u - u |v — | u " P u |v + | P u “ u |v ( 3 - 5 ) 

and 

| D t“ N - °t“|H i IV" - D t pN " u |H + l pND t“ - ‘V'Ih • ( 3 - 6 > 

A 

Lemma 3.1 and the fact that u has been assumed to be a strong solution 
imply that | P^u - ujy + 0 and j P^D^u - D^.uj^ + 0 as N + ® for each 
te[0, T]. 
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Since u satisfies (3.1) and (3.2) and since u satisfies (2.9) and (2.10) 

-w M 

then for 0 eV we have 

<M N Q D t 2 (u N - P N u), e N > H + aV - P N u, *0 N ) 

= < - P N u), e N > H + <(M 0 - mJ)D 2 u, 6 N > h 

+ a N (u-P N u, e N ) + a(u, e N ) - a N (u, e N ) 

+ b N (t)(u N - P N u, 0 N ) + b N (t)(P N u - u, 0 N ) 

+ b N (t)(u, 0 N ) - b(t)(u, e N ). 

, * , . v * M m m am a 

Let = (u^(t), u^j) = P u, let v = (v (t), v ) = u - and choose 

"N n *N W N T . 

0 = D^v eV . Then 


V\ + aV, D t ; N ) 

= <MqD^(u - P N u), D t v N > H + <(M q - M^u, D t v N > H 
+ a N (u - P N u, D^) + a(u, D^) - a N (u, D^v") 


+ b N (t)(v N , D t v N ) + b N (t)(P N u - u, D t v N ) 


+ b N (t)(u, o t ; N ) - b(t)(u, D t ; N ) 
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or 


1 n r^uNn A N n A N X . N, A N A N U 

• V >H * 3 (v * v » 


>mN/t D Ns n 2 A n A N K . WM M Nv n 2 A n A N. 

<Mq(I - P )D t u, D t v > H + <(Mq - H 0 )D t u, D t v > H 


+ D t a N ((I - P N )u, v N ) - a N ((I - P N )D t u, v N ) 

+ D t (a(u, v N ) - a N (u, v N )) - (a(D t u, v N ) - a N (D t u, v N )) 


, .N/,\/ a N pi A N \ , .N/.\/r)N A a n A N \ 
+ b (t)(v , D t v ) + b (t)(P u - u, D t v ) 


+ b N (t ) (u, D t v N ) - b(t)(u, D t v N ). 


Integrating both sides of the above expression from 0 to t, invoking 
hypotheses (HI) - (H3) on Q and using standard estimates we obtain 

min (trip m 2> m 3 )C | D t yN | H + | yN | 


✓ ft r , m 2 m 2 \ i / T d N x n 2 A i 2 , . n A Ni 2 

i/ 0 l max ( M l» M 3)|( I “ P ) D s u | H | D s v |H 


N \r.2 A , 2 


A Ni 2 


* max(|m - m N | 2 , |p - p"|')|D^u|^ + |D,/v"|‘ 


Ni2,i„2-.2 


‘Ni2 


+ M||(I - P N )D s u| 2 + 1^(2 + )E i _ EI-|-|D s u|- 


N.2. 


M^V + C 0 |a (q ^| 1 |^|^C 0 |^)| 1 |D s /|^ 


.No i A N i 2 


Ni 2 


+ C §|cr( q N )|f| (I - P")5|$ + |D s v"|j5 


N \ A i 2 


‘Ni 2 


+ C 1 |a(q) -a( q N )|2|u|2 + JD^} ds 
+ |U - pN ) u |v + C 2 1 | V + | ^ " pH ^\v 


N 1 2 1 A 1 2 


+ r- ei - Err u r, + c, v^r, + ei - ei 


V T 3 V 


‘Ni2 


,N,2, a .2 
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where Cq and Ci are constants which are independent of N and C2 and 
C3 may be chosen arbitrarily. Choosing C2 and C3 sufficiently 
small (C2 + C3 < min (mi, m2, m3)), we obtain 


I?! At) +/‘ 


where 


At) = K t {| (i - p N )i(t)|5 + | (i - P N )*|j; 

+ 1“ - n N |. 2 r»(t)|? ♦ |EI - Ei"|f|*l§ 

+ /S(|(I - pN ) D s“< s )|H + mx (| m - ” N | 2 . |p - P N |f) 

• | D s" u(s >|h + |» - P N )D s u(s}| 2 + | El - El N |f|D s S(s)|2 
+ |«(s, • ; q N )jf| (I - P N )'u(s)|^ + | 0 (s, • ; q) 

- o(s, • ; q N ) j f | u (s ) | y} ds 


A s ) " k 2^ k 3 + |°( s > * > q N ) | 

where K-j , i = 1, 2 , 3 are constants independent of N. Using the fact 
that + q in Q as N + « and u(q) is a strong solution together with 
Lemma 3 . 1 , the compactness of Q and Hypothesis (H 3 ) we have that 0 
and Kq is uniformly bounded in L2(0, T) as N + «. An application of 
the Gronwall inequality therefore yields 

In ^ M 1 2 1 A Ni 2 „ .. 

| D t v H + v V * 0 as N + ® 

for each te[0, T], Consequently |D t u N - D t P N u| H ^ 0 and |u N - P N u| v 0 
as N + » which together with ( 3 . 5 ) and ( 3 . 6 ) proves the theorem. 
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The continuous dependence results given in Theorems 3.2 and 3.3 
below can be verified using arguments similar to those used in the proof 
of Theorem 3.1. 

Theorem 3.2 : Let {q k } c 2 with q k q* as k If for each N fixed, 

u (q) denotes the solution to (3.1), (3.2) correspondingly to q then 

u N (q k ) + u N (q*) in V and D t u N (q k ) + D t u N (q*) in H for each te[0, T] as 

k “• That is, the mapping q + (u N (q), D t u N (q)) from Q into V x H is 
continuous for each N = 1, 2, ... 

Theorem 3.3: Let { q k } c 2 with q k * q* as k » and let u(q) denote the 

solution to (2.9), (2.10) corresponding to q. If u(q*) is a strong solu- 
tion then u(q k ) u(q*) in V and D^.u(q k ) D^.u(q*) in H for each te[0, T] 

A A 

as k -»■ «>. The mapping q -»• (u(q), D t u(q)) is continuous from Q into V x H 
in neighborhoods of those qeQ for which u(q) is a strong solution. 

We are now prepared to prove our first major convergence result. 

Theorem 3.4 : For each N = 1, 2, ... fixed, problem (IDN) has a solution 

_N _N, 

q . The sequence {q } admits a convergent subsequence {q } with 

- N k - 

q q e Q as k co. if u (q) is a strong solution to (2.9), (2.10) 
then q is a solution to Problem (ID). 

Proof 

Theorem 3.2 implies that the mapping q > J(q; u N (q)) from Q into R 
is continuous. This together with the fact that Q is compact yields the 

existence of a solution q^ to problem (IDN). The existence of a conver- 
- N k — W 

gent subsequence {q }c{q } also follows from the compactness of Q. If 
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- N k - 

q + q as k + « and u(q) is a strong solution to (2.9), (2.10) then 
Theorem 3.1 implies 

_ - _ _n. A .A 

J(q, u (q ) ) = lim J(q ; u k (q K )) 

k->-« 

A 

1 lim J(q; u K (q)) = J(q; u(q)) 

k+°° 

for all q e Q and the theorem is proven. 

Turning our attention next to the other infinite dimensional aspect 
of Problem (ID), we introduce a second level of finite dimensional approx- 
imation to effect a discretization of the admissible parameter space. 

For each M = 1,2, ... let the sets Q^c^be given by^rM(Q) where 
the mappings^ satisfy 

(PI) + 2 is continuous. 

M 

(P2) 3 (q) + q as M + » uniformly in q for all q e Q. 

We assume further that the sets have the property 

( P 3 ) For each qe(A u(q) is a strong solution. 

We define a doubly indexed sequence of approximating identification prob- 

1 ems by 

(IDNM) 

Find q = (m,EI,p,ag)eQ M which minimizes J(q; u^) subject to 

aM 

J(q; u ) subject to u being the solution to (3.1), (3.2) 
corresponding to q. 

Note that if the sets are of dimension < °° then the optimiza- 
tion in Problem (IDNM) is simply over a compact subset of the space 
R kM and is subject to finite dimensional constraints; a computa- 
tionally tractable problem. 
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The convergence arguments now go as follows. Property (PI) and Q 
compact imply that the sets Q M are compact as well. For each N = 1, 2, 
... and each M = 1, 2, ... Problem ( I DNM ) therefore has a solution 

q M . For each M = 1, 2, ... fixed, the sequence {qjj}^ admits a con- 

- N k - N k 

vergent subsequence {q M } with q M + q M as k -»• ~. Recalling property 
(P3) and arguing as we have in the proof of Theorem 3. 4, we conclude that 


q^ is a solution to the problem of minimizing J(q; u) over Q*^. Since 

Q M = ^(Q) there exists q M e Q such that cf M = ^ M (q M ). Now (q^cQ 

and Q compact imply the existence of a convergent subsequence {q^ } 

J 

- . M. - 

with q^ + q as j + «. Property (P2) then implies that^ J (q., ) -»■ q as 

J _ J* 

J + “ and hence that q M + q as j •»■ <» as wel 1 . Then 

j 


J(q M ; u) <_ J(q; u) 
j 



M. M. 

and Q J = & J (Q) yield 


_ M. A M . 

J ( q M. ; u ( q M.)) - J( q> ; u( ^ J (q)> q eQ * (3.7) 

J J 

Taking the limit as j + «> on both sides of (3.7) and recalling Theorem 
3.3 and Property (P3),we obtain 

J(q; u (q ) ) 1 J(q; u(q)) qeQ 
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from which we conclude that q is a solution to Problem (ID). In summary 
we have 

Theorem 3.5 : Each of the approximating identification problems 

( IDNM) has a solution . From the tableau {q|j} a sequence {q”(k ) } can 
be extracted with q*(k ) q", a solution to Problem (ID), as k + <=°. 

Typically an appropriate choice for the sets Q M are the spaces of 
linear or cubic interpolator spline functions with the mappings^ 
being constructed from the usual interpolation operators. If the ele- 
ments in Q are sufficiently regular, it is not difficult to verify that 
Properties (PI) - (P3) are satisfied for this approximation. A detailed 
discussion and several examples of this particular choice for the Q M in 
the context of inverse problems for parabolic systems can be found in 
[6]. Similar results for identification problems involving the estima- 
tion of functional parameters in beam equations with simple boundary con 
ditions (clamped, simply supported, cantilevered, etc.) are presented in 
[5] and [14]. 
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SECTION 4 


NUMERICAL RESULTS 


In this section we discuss numerical results for a variety of ex- 
amples. Although the analysis presented in the previous sections was 
based primarily on a simple example involving a cantilevered beam with 
tip mass, only minor modifications would be required so as to make our 
general approach applicable to a broad class of inverse problems for beam 
vibration. Some of these will be considered and outlined below. 

Although one of the major features of our scheme is its ability to 
identify or estimate functional (spatially and/or temporally varying) 
parameters, our numerical findings for this important class of problems 
are, at present, incomplete. For this reason, this report includes exam- 
ples involving the identification of constant parameters only. Our 
results for functional parameters will appear in a forthcoming paper. 

In all of the examples below, the fits were based upon artificially 
generated observations. By this we mean that so-called "true" values for 
the parameters were selected and a modal based Galerkin method was used 
to generate the solution to the resulting system of equations from which 
the sampled displacement measurements {u(t 1 -, x j)} were obtained * 

j=l,...,v 

Results for examples involving fits based upon actual experimental data 
will be discussed elsewhere. 
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The finite dimensional optimization problems were solved using an 
iterative steepest descent routine, ZXSSQ, from the IMSL Library. The 
algorithm used is the one suggested in the papers by Levenberg [15] and 
Marquardt [18], The finite dimensional second order initial value prob- 
lem given by (3.3), (3.4) which has to be solved at each iteration to 
compute the value of the least squares payoff functional J and approxima- 
tions to its gradient and the corresponding Jacobian matrix is integrated 
using a variable order Adams predictor corrector method (IMSL routine 

DGEAR). The system (3.3), (3.4) did not in general demonstrate stiff 
behavior. 

The integral inner products in the generalized mass and stiffness 
matrices and the generalized Fourier coefficient vectors for the external 
loads and initial data were computed using a composite two point Gauss- 
Legendre quadrature rule. The four subinterval support of the cubic 
B-splines and their derivatives leads to 7-banded matrices and conse- 
quently contributes to the overall computational efficiency of the 
method. 

All examples were run on the IBM 3081 processor at Draper Labora- 
tory. 

Example 4.1 

In this example we consider a cantilevered beam of length 1.0 with 
an attached tip mass. We seek to identify the spatially invariant stiff- 
ness El and linear mass density p of the beam and the magnitude m of the 
tip mass. We assume that the system is initially at rest ($ = 0, ^ = 0) 
and then excited by the distirbuted lateral load along the beam given by 

f(t, x) = e x sin 2-rrt 
and the point force applied to the tip mass 

g(t) = 2e' 2t . 

We assume that there is no axial loading, or, a = 0. 
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Displacement observations were generated by taking El = 1.0, p = 3.0 
and m = 1.5 to be the true values of the parameters. Measurements were 
taken at positions xj = 0.25j, j = 2, 3, 4, at times tj = 0.5i, i = 

1, 2, ... 10, from a solution to the system (2.1) - (2.5) generated using 
a Galerkin method and the first two natural modes of the unforced system 
(see [25]). The "start up" values for the steepest descent routine were 
EIq = 0.7, po = 2.7 and mo = 0.7. The final converged values for 

M M — |S| 

El , p , and m together with the residual sum of squares J and re- 
quired CPU times for various values of N are given in Table 4.1 below. 

Table 4.1. 


N 

if 

-N 

P 

-N 

m 

3* 

CPU 

(m:s) 

2 

1.00057 

3.04455 

1.48957 

0.12 x 10" 3 

0:9.19 

3 

1.00067 

3.01256 

1.49707 

0.11 x 10" 3 

0:22.10 

4 

1.00027 

3.00922 1 

1 

1.49721 

0.11 x 10" 3 

0:57.60 

5 

1.00016 

I 

2.98936 

1.50262 

0.11 x 10" 3 

1:22.52 

6 

0.999T2 

i 

2.99720 

l 

1.49952 

0.11 x 10' 3 

1 

2:52.76 


Example 4.2 

We consider the system described in the previous example. We assume 
that it is initially at rest and then excited by the distributed load 

f(t, x) = e X sin 2ut 
along the beam and the point force 

g(t) = 2e _t 


33 




acting on the tip mass. We also assume that the entire system is sub- 
jected to a base acceleration which is given by 


a 0 (t) 


( 1 0 < t £ 1.5 

( 0 otherwise. 


(4.1) 


The internal tension resulting from the axial load (see [26], [27]) is 
given by 


o(t, x) = -a Q (t)(pU - x) + m). 


We are interested once again in estimating the stiffness El, linear mass 
density p,and the magnitude of the tip mass m. The true values of the 
parameters were taken to be El = 1.0, p = 3.0, and m = 1.5 with the 
reference solution being generated using the first two natural modes of 
vibration for the unforced, unaccelerated system. Displacement measure- 
ments were taken at positions xj = 0.75, 0.875, 1.0, j = 1, 2, 3, at 
times t-j = 0.5i , i = 1, 2, ..., 10. The start up values for the itera- 
tive search routine were taken as EI 0 = 0.7, po, and mo = 1.7. Our 
results are summarized in Table 4.2. 


Table 4.2. 


N 

ET* 

-N 

P 

-N 

m 

j* 

CPU 

(m:s) 

2 

1.00057 

3.09966 

1.47928 

0.17 x 10" 4 

0:18.01 

3 

1.00121 

3.06360 

1.48727 

0.18 x 10" 4 

0:35.93 

4 

1.00092 

3.04144 

1.49207 

0.19 x 10" 4 

1:26.19 

5 

1.00057 

3.03063 

1.49413 

0.19 x 10 -4 

4:15.87 

6 

1.00117 

3.03186 

1.49436 

0.35 x 10" 4 

5:21.34 
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We note that although strictly speaking the convergence theory 
developed in Section 3 requires that aoeC 1 , the scheme performed 
satisfactorily with ao given by (4.1) above. 

Example 4.3 

In this example we consider a free-free beam of length 1.0 with an 
attached tip body at each end (see Figure 4.1). 



Figure 4.1. 

The tip bodies are assumed to have known mass properties which are given 
by: 

Tip Body 0 (at x = 0): 

mo = 0.75, co = 0.1» So - ir/ 6 , Jo = 0.6 
Tip Body 1 (at x = 1): 

mi = 1.5, ci = 0.2, Si = tt/3, Ji = 0.4 
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where for tip body i, i = 0, 1, itk , c.. , 6., and J. are respectively its 
mass, the distance from its center of mass to the tip of the beam, its 
mass center offset as measured from the extension of the longitudinal 
axis of the beam, and its moment of inertia about its center of mass. 

The equations describing the transverse displacement of the beam and 
the translational and rotational equilibrium of the tip bodies are given 
by (see [25], [27]) 


pD t u + D x EID x u = f ’ 


xe (0, 1), te (0,T) 


(4.2) 


rngD^u - hIqCq cos 6qD^D x u + D x EID x u = 0» x = 0, te (0, t) 


2 2 2 2 
-m 0 c 0 cos 6 qD^.u + (J Q + m 0 c 0 )D t D x u - EID x u 


= 0, x = 0, te (0, T) 


m^D^u + m-^Cj cos <$jD^D x u “ ^ x E ^x u = x = 1, te (0, T) 


-Hn^c ^ cos 5 i^t u + ^1 + m l c l^t^x u + EID x u 


= 0, x = 1, te (0, T). 


The initial conditions are of the form 


u = <f> 


D t u = * 


xe (0, 1), t = 0 


(4.3) 


Letting H = x H®(0, 1) with inner product 


<(n» <|»), (?, ^)> H = n T ^ + <<f», ^>g 


and 

V = {}eH:J = ( n , +), <f»eH 2 (0, 1), n = (<f>(0), D*(0), *(1), D<j>(1)) T } 


36 



with inner product 


<(J) , <(j) , Tp> £ • 

The weak form of the system given by (4.2), (4.3) above becomes 

<M Q D 2 u(t), e> H + a (u (t ) , 0 ) = <F 0 (t), 0 > H te (0, T), 0 e V (4.4) 

u(0) = l D t u(0) = i (4.5) 

where 

u(t) = ( (u (t, 0 ), D x u(t, 0), u(t, 1), D x u(t, 1)) T , u(t, • ))eV , 

Mg(n * <f>) = (Mgn, p<}>) 

with 



a: V x V + R given by 

a(<f>, ip) = <E ID 2 (f> , D 2 i|)> 0 , 

F 0 (t) = ( 0 , f(t, .)) 

J = (4(0), D<j>(0), 4,(1), D<j>(l)) T , 4 ) 
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and 


♦ = ((*( 0 ), MO), 4.(1), MD) T , 4 ,). 

The appproximation scheme for the system (4.4), (4.5) is then constructed 
in essentially the same manner as it was in Section 3. 

In this example we seek to identify El and p. The system was as- 
sumed to be initially at rest and then excited via the input disturbance 

f(t, x) = 10e X sin 2irt. 

Observations at positions Xj = 0.25j, j =0, 1, ..., 4, at times tj = 
0.2i, i = 1, 2, ..., 10, were generated using the first six natural modes 
of the system; four flexible modes plus rigid body rotation and transla- 
tion. The true values of the parameters were assumed to be El = 1.0 and 
p = 3.0 with start up values taken as EIq = 0.7 and po = 2.7. Our 
results are given in Table 4.3. 


Table 4.3. 


N 

if 

-N 

P 

J* 

2 

0.99567 

3.00092 

O 

• 

no 

0 

X 

i—* 

0 

1 

4* 

3 

0.99374 

2.99900 

0 

• 

h- * 

CT> 

X 

l— » 

0 

1 

4 

0.99849 

2.99798 

LO 

O 

r— 1 

X 

CO 

r^. 

• 

0 

5 

0.99888 

2.99910 

0.28 x 10" 5 
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Example 4.4 

In this example we estimate the flexural stiffness El and linear 
mass density p for a cantilevered beam of length 1.0 with an attached tip 
body (see Figure 4.2). 



Figure 4.2. 

We assume that the system is initially at rest and then acted upon by the 
distributed load 


f(t, x) 

and base acceleration 


20e" 2t e " 2 °( 1-x ) 


a 0 (t) 


1.0 


0 

1.0 

0 


0 < t < 1.5 
1.5 < t < 3.0 

3.0 <_ t < 4.0 

4.0 < t. 


39 


The equations of motion are given by (see [26], [27]) 


pD^u + EID^u = “ a q D x ( p ( 1 - x) + m)D x u + f xe(0,l), te(0,T) 

mD t u + me cos 6D^D x u - EID^u = ma^u, x = 1, te(0,T) (4.6) 

me cos 6D 2 u + (J + mc 2 )D^D x u + EID 2 u = me cos 6a 0 D x u + me sin Sa Q , 

x = 1, te(0,T) 

with boundary and initial conditions 

u = 0 D x u = 0 x = 0, te[0, T] (4.7) 

and 


u “ <f> - ip xe[0, 1], t = 0 (4.8) 

respectively where the quantities m, c, 6, and J are as they were defined 
in Example 4.3. Once again only minor modifications (similar to those 
outlined in the previous Example) are necessary so as to make the theory 
presented in Sections 2 and 3 applicable to the system given by (4.6) - 
(4.8) above. Taking the true values of El and p to be 1.0 and 2.0 re- 
spectively and setting m = 4.0, c = 0.2, 6 = tt/ 3, and J = 0.4, displace- 
ment observations at positions xj = 0.75, 0.875, 1.0, j = 1, 2, 3, at 
times t-j = 4.0, 4.5, 5.0, i = 1, 2, 3, were generated using the first 
three natural modes of the unforced, unaccelerated system. The start up 

values were taken as EI 0 = 0.7 and po = 2.5. The final converged 

— N “ N 3 m 

values for El and p together with the residual sum of squares iT for 

several values of N are given in Table 4.4 below. 
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Table 4.4. 


1 

FT* 

-N 

P 


1 

1.00046 

2.07731 

0.4 x 10‘ 3 

9 

1.00135 

2.06652 

0.22 x 10" 3 

■ 

1.00117 

2.05180 

0.23 x 10" 3 

1 

1.00086 

2.04776 

0.26 x 10" 3 
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